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1. Introduction 



In this proceedings contribution, we summarize the computation of the fragmentation of a 
gluon into a spin-triplet 5-wave quarkonium via the ^S^l^ channel in relative order v^, where v is 
the velocity of the heavy quark Q or heavy antiquark Q in the quarkonium rest frame. We refer the 
reader to Ref. [|I|] for details of this calculation. The expression 2^+1/ W the standard spectroscopic 
notation for the spin s, orbital angular momentum I, total angular momentum j, and color c = 1 
(singlet) or 8 (octet) of the QQ pair that is created at short distances and that evolves into the heavy 
quarkonium. 

Our computation make use of the nonrelativistic QCD (NRQCD) factorization approach [^. 
That is, we express the fragmentation functions as sums of products of NRQCD long-distance ma- 
trix elements (LDMEs) and short-distance coefficients. The focus of our calculation is the short- 
distance coefficient that appears in the order v"^ contribution to the fragmentation function for the 
^5'['' channel. The short-distance coefficients for gluon fragmentation in this channel in relative 
orders v'' and have already been computed in Refs. ^ and [Q], respectively. The contribution 
of relative order v'^ is interesting theoretically because it is at this order that the ^sf^ fragmentation 
channel first develops soft divergences. These soft divergences are ultimately absorbed into the 
NRQCD LDME for the evolution of a ^S^f^ or ^pf^ QQ pair into a spin-triplet ^-wave quarkonium 
state. Consequently, the ^sf^ contribution in order v"* and the ^S^f^ and ^pf^ contributions at the 
leading nontrivial order in v are related by logarithms of the factorization scale. Since the ^S^f^ con- 
tribution to the fragmentation of a gluon into a is known to be significant phenomenologically, 
it is important to compute the ^sf^ contribution in order v^, as well. 

An important technical issue in this calculation is the appearance of both single and double 
soft poles in the dimensional-regularization parameter e. This is the first NRQCD factorization 
calculation in which double soft poles have appeared. In order to effect the matching between 
NRQCD and full QCD that removes these poles, it is necessary to work out, for the first time, 
two-loop corrections to NRQCD operators and the corresponding ultraviolet renormalizations. 

We employ the CoUins-Soper definition of the fragmentation function for a gluon fragment- 
ing into a quarkonium. We compute the fuU-QCD fragmentation functions for a gluon fragmenting 
into free QQ states with various quantum numbers. We then determine the NRQCD short-distance 
coefficients by comparing the fuU-QCD fragmentation functions with the corresponding NRQCD 
expressions, making use of the NRQCD factorization formulas. 

The remainder of this paper is organized as follows. In Sec. ||, we present the NRQCD factor- 
ization formulas through relative order for the fragmentation functions for a gluon fragmenting 
into aJ/Y- The results of the fuU-QCD calculation for the corresponding fragmentation functions 
for free QQ states are given in Sec. |[ In Sec. ^, we compute the relevant NRQCD LDMEs for free 
QQ states analytically in dimensional regularization, and we determine the evolution equations 
for the LDMEs. In Sec. |5|, we compute the short-distance coefficients by matching the NRQCD 
and fuU-QCD results. Sec. |6| contains estimates of the relative sizes of the various fragmentation 
contributions to the J/y hadroproduction cross section. We give a summary of our results in Sec. ^. 
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2. Factorization formulas 



We denote by D[g — )• 7/i/a](z,/Ia) the fragmentation function for a gluon fragmenting into 
Si J/Y- Here, jU.^ is the factorization scale, z = P^/k^, and P and k are the momenta of the 
J/Y and the fragmenting gluon, respectively. We define light-cone coordinates for a four- vector 
V = {V+,V-,V±) by = {V^±V^)/V2. We make use of the Collins-Soper definition of 
D[g -^J/\ff] {z, IJ-a) which is gauge invariant, and we regularize soft and ultraviolet divergences 
dimensionally, taking d = 4 — 2e space-time dimensions. 

The NRQCD factorization formula for D[g — ^ J/ i//](z,/Xa) is given by 



Z)[g^7/VA](z) = £c/„(z)(0|^/'^|0), 



(2.1) 



where the (0|^,f/'^|0) are NRQCD LDMEs, the J„(z) are short-distance coefficients, and we have 
suppressed the /^a dependences in <i„(z) and {0\i?i^'^\0). Since the d„{z) are independent of the 
hadronic final state, we can write the fragmentation function for a gluon fragmenting into a free 
QQ state as 

D[g ^ !2G](z) = £4(z)(0|^2e|o). (2.2) 



By computing the left side of Eq. (2.2) in full QCD and comparing it with the right side, we can 
determine the (z). In the remainder of this paper, we denote the order- v*^ contribution to D[g — )• H] 
by Dk[g — )• H], where H can be either a quarkonium or a free- 22 final state. 
The operator LDMEs that we use in this paper are 



(0|< 



0) 
0) 
0) 
0) 

0) 

0) 

0) 
0) 



^(0|z^a'(-^O)V=^//V^Vz + H. c.|0), 



^{0\x''a'{-^D)'Y V^o'x + ^. c.|0), 



:Dfx\0), 



1 



JO\X^o'y V^<y'{D-g,,E + g,E-D)x 

-X^a'{D-g,E + g,E-D)Y ^H¥''(y'x\0), 
{0\x^a'T"Y^H¥^o'T''x\0), 

{o\xH-^DYa"r'Y vH-^DYo''r'x\o), 



(2.3) 
(2.4) 
(2.5) 
(2.6) 

(2.7) 

(2.8) 
(2.9) 

(2.10) 



where ^.v = \/47ra^., and x are two-component (Pauli) fields that create a heavy quark and a 
heavy antiquark, D is the gauge-covariant derivative, and E is the chromoelectric field operator. 
^H(P) = Lx l^(^) +X){H{P) +X\ is a projection onto a state consisting of a quarkonium H, 
with four-momentum P, plus anything. ^h{p) contains a sum over any quarkonium polarization 
quantum numbers that are not specified explicitly. The symmetric traceless product is defined 
by A^'B}^ = ^{A'Bj + A^B') - ^ S'^A'^B'', and the antisymmetric product is defined by A^'B^^ = 
^{A'BJ -AJ'B'). 
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We now give the NRQCD factorization formulas for gluon fragmentation into a through 
relative order v^. The contributions to D[g — )• in relative orders v°, v^, and are 

Do[g^J/w\ = ^ofe^G<2('5W)](0|<'"^(35W)|0), (2.11) 
D2[g^j/W\ = ^2fe^G<2('5W)](0K'"^(^5W)|0), (2.12) 
DAg^J/w] = d,[g^QQ{'sf)mir',l''{'sfm. (2.13) 

The short-distance coefficients d^lg QQCS^^^)] and diig QQCS^^^)] were calculated previ- 
ously in Refs. (jsj, ^ and Ref. respectively. The short-distance coefficient do[g — QQ{^Sq^)] 
is related by an overall color factor to do[g — QQ^Sq^)], which was calculated in Ref. [^. The 
contribution to — in relative order is 

D4[g^j/v] = {d4Ag^QQCs?)]+d4A8^QQCs?)]}{o\^J''es[^)\o) 

+ do[g^QQep^'^)]m'j''ep^'^m 

+ do[g ^ QQesP)m^/''Csf^)\0). (2.14) 

The short-distance coefficient do[g — )• QQ^sf^ )] was calculated previously in Refs. [||, ^ ||, ^]. The 
short-distance coefficient do [g — QQC'P^^^ )] is related by an overall color factor to the sum over J of 
the short-distance coefficients do[g — QQCP^^)] that were calculated in Ref. [^. The computation 
of the combination of short-distance coefficients d^^g QQCS^I^)] +d4Ag — )• QQC'S^^^)] is the 
main goal of this work. ^ 



3. FuU-QCD calculations 



We compute the fragmentation functions for free QQ states in full QCD in J = 4 — 2£ dimen- 
sions, multiplying the square of the dimensional-regularization scale /i by the factor e^^/ (47r) that 
is appropriate to the modified-minimal-subtraction (MS) scheme. Here, is the Euler-Mascheroni 
constant. 

The results for Dolg QQ{^sf^)] and D2[g QQ{^P^^^)] at leading order (LO) in a, and v 



are 



Do[g 
Diig- 



\47l 



8«2g2 

(J-l)m5 4Nc 
1 



5(1 
4 



z), 

l-£)r(l + e) 



(3.1) 



471 



4m2 



2£iR 



5(l-z)+/(z) 



(3.2) 



'in fact, the color-singlet contribution to D4[g J/w] is a linear combination of the LDMEs (0|^;'''("'5['')|0), 
for ( = 1-3. We express the contribution of {0\ff^^^(^s'^^^)\0) in terms of the other two LDMEs by making use of the 
NRQCD equations of motion, and we use the vacuum-saturation approximation to replace the {0\ff'^{'^(^s'^i^)\0) with a 
single LDME, (0|^'''''(^5['')|0), making an error of relative order ||l^. 
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where Nc is the number of colors, q is half the relative momentum of the Q and Q in the QQ rest 
frame, and the finite function f{z) is defined in Eq. (5.20) of Ref.|[l]]. The subscript "IR" in e^^ 
indicates that the pole is infrared in origin. 



The result for D^[g QQQS^I"^)] at LO in a, is 



■D,[g^ eeCs™ +/[5i2] +i[S2] 



where /[S^], ^[^i], and /[^i] contain the double- and single-pole contributions: 



^ 5(1 -z) ^ 



2eiR 



1 



2 m2 



(l_^)l+4e 



+ 



l-z 



l+2e 



2£ir(1 



U+4e 



/ 8«, \ 

\3nm^) 16A^2 {^d - 1)171^ \Ati 



2 \ e 4 



d-l 



4m^ 



I[Si 



nsi] 

( 8a, 



5_ 

2£iR 

2 ^t2 



4m2 

2 N e 



2e 



,-2+2e/ 



(l_^)-4er2(l+£) 

48(1 -£) 
■ + Oie). 



(3.3) 



) l!(i±^E!|i_M(,_,),6-2.-.^--2e^), (3.4, 



(3.5) 



) \6N^ \4k J {d-l)m^d-\ 
Here, Ti and To are functions of z that are defined in Eq. (5.29) of Ref.[|T]|. D^lg QQ{^ S^l^)f'^'^'' 



in Eq. (3.3) is a finite contribution, which we evaluate numerically in J = 4 dimensions. 



4. NRQCDLDMEs 



_ T Ml 

In this section we tabulate our results for the NRQCD LDMEs for free QQ{^S\^ ) states that 
are relevant through relative order v^. 

The free QQ matrix elements at order af are normalized as 

(4.1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 

where there is an implied sum over final-state polarizations, and the superscript {k) indicates the 
order in a,. 

In order a, and order a^, the relevant LDMEs mix: 

-l\A^2_i 1 

'MS Q^,^2 \ oct: ) A All 1 ,^^\^ir ' r-i 



(0|. 




)|0)(°) 


= (A^.'-i), 




(0|. 




)|0)(°) 


= 2(J-l)A^„ 




(0|. 




)|0>(°) 


= {d-l)iN^- 


1), 


(0|. 




)|0)(«) 


= g'mf^''' 


)(^5S"1)|0)«'), 


(0|. 




)\or 




^(^5S"1)|0)(0), 






)|0)(°) 


= qHd-l){N^ 


-1), 



(UFo ^'"^MS - 3;j,„2 U&rJ 4A^2 j_i^"I^4.i 



(0) 



(4.7) 
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(0|^?««""'>(^M")|0)™ 



(oi^f^i"Vsri)io)£ 



)|0> 



(4.8) 



(0) 



(4.9) 



Here, the subscript MS indicates that we have removed the ultraviolet poles in e by using the 
minimal-subtraction procedure, with the choice of scale that is appropriate to MS subtraction. 

The renormalization subtractions for (0|^^^' ^(^5f^)|0)'^' involve both Bom diagrams for the 
double-pole counterterm and one-loop diagrams for the single-pole counterterm. In the minimal 
subtraction procedure, one subtracts only pure pole contributions in the ultraviolet-divergent sub- 
diagrams. It is essential for the consistency of the minimal-subtraction program to treat factors 
that are external to those divergent subdiagrams, such as angular-momentum projections involving 
external momenta, exactly in J = 4 — 2£ dimensions. A failure to follow this procedure in the 
presence of poles in J = 4 — 2e of order two or higher can result in the appearance of uncanceled 
poles in e in the NRQCD short-distance coefficients. We refer the reader to Ref. [[T]] for details. 

The dimensional-regularization scale ^ can be identified with the NRQCD factorization scale 
/Ia. It then follows from (3fai/(ilog(/iA) = —leas + 0{a'^) that the renormalization-group evolution 
equations for the LDMEs are 





JlogjUA 

d 

"5? log Ma 
d 

log Ma 



^0|^eeeM")^3p[8] 



(1^ 

MS 



8 a, 



3Tinfi AN} d - 1 



(0|, 



yeGe5W)^3^[l])|Q^(0) 



ml 



,ee('si") 



ill 

MS 

{2)_ 
MS 



8 a, Nt 



?>nnP- 4Nc 
8a, N}-4 



t^o|^eeepi«i)(3p[8])|o^(o)^ 



(4.10) 
(4.11) 
(4.12) 



durrfi 4Nc 

The result in Eq. ( 4.10| ) agrees with the corresponding results of Refs. [11, |l^]. Equations ( |4.1 1| ) 
and ( 4.12 ) agree with the result in Eq. (B19b) of Ref. at the leading nontrivial order in v and 
with the corresponding result in Ref [12], but disagree with the corresponding result in Ref. [fTJ. 



5. Short-distance coefficients 



By making use of the results of Sees. ^ and ^ and the fvee-QQ versions of the NRQCD factor- 
ization (matching) equations ( 2.11 ), ( 2.12 ), and ( p.l4[ ) in Sec. |^ we obtain 



^4 fe^eeC^w + 

-xlog — + log^— + 
3 2m 2m J 



3n{d-\yN^m'^ 
1 



log(l 



104-29z-10z2 7[z + (l+z)log(l-z) 



24 

+ i±^(31-6z- 



■ + 
36 



+ ■ 



1-z 
l-2z)(8-5z) 



2z2 

log(l-z)-^^39-6z+^_^ 



log 



+ 

Ma 
2m 



logz 
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Figure 1: The color-singlet short-distance coefficients t/ofe -> QQ{^S^l^)]^^\z), d2[g QQ{^ S^l^yp'\z) and 



d^lg — > QQ(j'S^l^)]^"^^^{z), which are defined in Eqs. ( 2.11 ), ( 2.12 ) and (yA), respectively, as functions of z- 
The scaling factors are (^,^,^) = (10^^ x a?Jnr',\Q-^ x a'^Jnr', 10"^ x a?Jm''). 



+ ■ 



13 -7z 



l0£ 



:^-l0gg)l0g(l-z)-Li2(z) 



Here, d^[g QQCS^^)]^""^^ = D^g QQCS^I^)]^''''^ {0\fff^^'^^'^\^s[^)\0)^°^ is evaluated by car- 
rying out the integrations over the phase space numerically. The results of the numerical integration 
are shown in Fig. [I| One can find details of the matching procedure in Sec. 7 of Ref. 



+ 0(8) 



(5.1) 



6. Relative sizes of the fragmentation contributions 

We now make estimates of the relative sizes of the contributions of the various fragmentation 
functions to the cross section for //i// production in hadron-hadron collisions. At large pr, the 
differential cross section for a gluon can be approximated as dOg/dpj- <x \/pj, where is a fixed 
power. Then, the fragmentation contribution to the production rate becomes dOjJ^/dpj oc 

Ik{D), where Ik{D) = J^dzz'^ D{z) ijl]]. Hence, we can obtain a rough estimate of the relative 
contribution of a fragmentation process to the cross section by computing Ik{D). At large p^, the 
cross section at LO in is dominated by fragmentation of a gluon into a. J/y with z = 1 [p^]. 
The next-to-leading order (NLO) k factor for this channel is essentially independent of p^ for the 



cross sections that were measured at the Tevatron [|14|. Making use of these facts and the result 
for the ^sf' contribution to da/dpj at NLO in that appears in Fig. 1(c) of Ref. [15], we have 
determined that K" « 5.2. 



We can compare the relative contributions to dOjj^ of each of the three QQ channels that 
contribute to Z)[g — in order [Eq. (2. 14)]. For each channel, we compute h.2{D). Then, we 
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Table 1: Relative contributions to dOjJ^ in order v^. The first and second rows give h.2{d) times the LDMEs 
that are given in the text. We take m = mc = 1.5 GeV. For compatibiUty with Ref. [p^, we take a., = ttj ('«r)> 
where nij = J pj+4m^. We choose the point pj- = 20 GeV, which impHes that as{mj) — 0.154. 



/^•((i) \ channel 












h.2{d)\nj^=m X 


lO^LDME 


1.65 


18.6 


32.7 


0.192 


h.2{d)\^i^=2m X 


lO^LDME 


1.65 


18.6 


43.6 


0.455 



multiply by the following LDMEs: (0|^/"^(^5['1)|0) = ni^{vY x 1.32 GeV^ {0\^J'^ {^P^^^ 
-0.109 GeV5, (0|^/"^(35fl)|0) =3.12 X 10-3 GeV^, {0\^J'^ {^sf)\0) =4.50x 10"^ GtW\ We 
have obtained these LDMEs from Ref. [[151] . In the case of the first LDME, we have used the 
generalized Gremm-Kapustin relation ^ (0|^/"^(35W)|0) =m4(v2)2(0|^/"^(35['l)|0) to obtain 
the matrix element of order v"^ from the matrix element of order v° that appears in Ref. |JT5|]. We use 
the value of rr?{v^) from Table I of Ref. ||l^]: rr?{\^) = 0.437 GeV^ at m = 1.5 GeV. The results 
of this computation are shown in Table |I[ 

We see from the second row of Table [I] that the QQi^S^l^) channel makes a small contribution 
to Dn[g ^ J /'[^f] at pj = 20 GeV, confirming that this channel is not important phenomenologically 
at the current level of precision. We also see that the QQi^sf^) and QQi^P^^"^) channels give 
comparable contributions to D^[g ^ J /\if\ at pj = 20 GeV. We estimate that the fragmentation 
contribution to do /dpj x B{J/y — > MM) at pj- = 20 GeV from the QQij'P^^^) channel is about 
6 X lO^^* nb/GeV [jl|], which is comparable to (about a factor of 2 larger than) the total NLO 
contribution in the QQ{^P^^^) channel in Fig. 1(c) of Ref. [15]. A more precise calculation of 
the fragmentation contribution in the QQ{^P^^^) channel will be necessary in order to determine 
whether it is the dominant contribution in that channel at NLO in at high p^ . 

7. Summary 

We have computed NRQCD short-distance coefficients for gluon fragmentation into a J/y 



through relative order v^. Our main result is the expression in Eq. (5.1 1 for (^41 [g — QQ^s'^^)]^^'^ + 



d4,2[8-^QQCS^i^)] ' which is the sum of short-distance coefficients for the ^S^l^ NRQCD LDMEs 
of relative order v"*. This is the first time that double soft divergences, as well as single soft diver- 
gences, have appeared in intermediate steps in the calculation of NRQCD short-distance coef- 
ficients. We have also computed the free QQ NRQCD LDMEs that appear when one uses the 
matching equations between full QCD and NRQCD to compute the NRQCD short-distance coef- 
ficients. Our perturbative calculations of these LDMEs involve both one-loop and two-loop cor- 
rections. The ultraviolet renormalization-subtractions for these LDMEs include Born and one-loop 
single-pole counterterm contributions and Born double-pole counterterm contributions. We have 
worked out the renormalization-group evolution of the renormalized LDMEs and confirm some 



previous results in Refs. [0, 1 1 , 12|, but disagree with a result in Ref. [11]. 



We have estimated the relative sizes of the contributions of the various channels to gluon 
fragmentation into slJ/y through relative order v^. The contribution in order of the ^S^^^ channel 



8 



Order-v^ Relativistic Corrections to Gluon Fragmentation into ^Si Quarkonium 



Jungil Lee 



to the cross section at = 20 GeV is about a factor of 2 (4) larger than the contribution in order 
when /Xa = m^. (Inic). In spite of this large enhancement of the fragmentation contribution in 
order v^, the corresponding contribution to the 7/ production cross section at the Tevatron or the 
LHC is not important at the current level of precision. 

In the ^sf^ channel, most of the contribution to the J/y production cross section at hadron- 



hadron colliders at large /j^ arises from the fragmentation contribution. This is true at LO in [13] 



and atNLO in [[l4|]. Our estimate of the fragmentation contribution to the ^Pj channel indicates 
that it is an important part of the contribution in that channel to the high-;?^ 7/ y/ production cross 
section at hadron-hadron colliders at NLO in a^. However, it will be necessary to carry out a more 
precise calculation of the fragmentation contribution in the channel in order to see whether it 
is actually dominant at high p^. . 
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